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Although meta-analyses of genome-wide association studies have 
identified >60 single nucleotide polymorphisms (SNPs) associ- 
ated with type 2 diabetes and/or glycemic traits, there is little 
information on whether these variants also affect a-cell func- 
tion. The aim of the current study was to evaluate the effects 
of glycemia-associated genetic loci on islet function in vivo and 
in vitro. We studied 43 SNPs in 4,654 normoglycemic participants 
from the Finnish population-based Prevalence, Prediction, and 
Prevention of Diabetes-Botnia (PPP-Botnia) Study. Islet function 
was assessed, in vivo, by measuring insulin and glucagon concen- 
trations during oral glucose tolerance test, and, in vitro, by mea- 
suring glucose-stimulated insulin and glucagon secretion from 
human pancreatic islets. Carriers of risk variants in BCLllA, 
HHEX, ZBED3, HNFIA, IGFl, and N0TCH2 showed elevated 
whereas those in CRY2, IGF2BP2, TSPAN8, saidKCNJll showed 
decreased fasting and/or 2-h glucagon concentrations in vivo. 
Variants in BCLllA, TSPAN8, and N0TCH2 affected glucagon 
secretion both in vivo and in vitro. The MTNRIB variant was 
a clear outlier in the relationship analysis between insulin secre- 
tion and action, as well as between insulin, glucose, and gluca- 
gon. Many of the genetic variants shown to be associated with 
type 2 diabetes or glycemic traits also exert pleiotropic in vivo 
and in vitro effects on islet function. Diabetes 62:2978-2983, 
2013 




In the last few years, genome-wide association stud- 
ies have substantially increased the knowledge of 
genetic variants predisposing to type 2 diabetes. Al- 
though the majority of these single nucleotide poly- 
morphisms (SNPs) seem to influence insulin secretion, 
few, if any, studies have assessed their simultaneous ef- 
fects on (3- and a-cell function in vivo and in vitro. The 
aim of the current study was to provide a comprehensive 
evaluation of the effects of genetic loci associated with 
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type 2 diabetes (1) and/or glucose and insulin levels (2) 
on islet function in vivo, in a large well-characterized 
population-based study from the western coast of Finland 
(the Prevalence, Prediction, and Prevention of Diabetes- 
Botnia [PPP-Botnia] Study), and in vitro, in human pan- 
creatic islets. Islet function was assessed by measuring 
insuUn and glucagon concentrations during an oral glucose 
tolerance test (OGTT). 



RESEARCH DESIGN AND METHODS 

Study population. The PPP-Botnia Study is a population-based study from the 
Botnia region of western Finland. Nondiabetic subjects (n - 4,654; fasting 
plasma glucose <7.0 mmol/L and 2-h plasma glucose <11.1 mmoI/L) were 
included in the current study (3) (Supplementary Table 1). 
Measurements. Blood samples were drawn at 0, 30, and 120 min of the OGTT. 
Insulin was measured using a tiuoroimmimometric assay (AutoDelfia; Perkin- 
Ehner, Waltham, MA) and serum glucagon using radioimmunoassay (MiUipore, 
St. Charles, MO). Insulin sensitivity index (ISI) was calculated as 10,000/ 
VCfasting P-glucose X fasting P-insulin X mean OGTTgiucose X mean 
OGTTinsuiin) (4). lusuliu secretion was assessed as corrected insulin re- 
sponse (CIR) during OGTT {CIR = 100 X insulin at 30 min/[glucose at 30 minX 
(glucose 30 min - 3.89)]} (5) or as disposition index (DI), i.e., insulin secretion 
a4justed for insulin sensitivity (DI = CIR X ISI). 

Genotyping. Genotyping was performed either by matrix-assisted laser 
desorption-ionization time-of-tiight mass spectrometry on the MassARRAY 
platform (Sequenom, San Diego, CA) or by an allelic discrimination method 
with a TaqMan assay on the ABI 7900 platform (Applied Biosystems, Foster 
City, CA). We obtained an average genotyping success rate of 98.1%, and the 
average concordance rate, based on 10,578 (6.5%) duplicate comparisons, was 
99.9%. Hardy-Weinberg equilibrium was fiilfllled for all studied variants (P > 0.05). 
Glucose-stimulated insulin and glucagon secretion in human pancreatic 
islets. Glucose-stimulated insulin and glucagon secretion were performed in 
high (16.7 mmol/L) and low (1.0 mmol/L) glucose concentration in the medium 
as described previously (6) and were available from 56 nondiabetic human 
pancreatic islet donors. Given the limited number of donors for genetic 
analyses, the relationships between absolute genotypic means of insulin and 
glucagon concentrations and the 95% CI were plotted separately for different 
genotype carriers rather than contrasting mean differences in phenotypes 
between genotypes. 

Statistical analysis. Variables are presented as mean ± SD and, if not nor- 
mally distributed, as median (IQR). Genotype-phenotype correlations were 
studied using linear regression analyses adjusted for age, sex, and BMI. 
Nonnormally distributed variables were logarithmically (natural) transformed 
for analyses. All statistical analyses were performed with IBM SPSS Statistics 
version 19.0 (IBM, Armonk, NY) and PUNK version 1.07 (7,8). 

RESULTS 

Effect of SNPs on metabolic variables 

Insulin secretion. We observed the strongest effect on 

reduced insulin secretion for the MTNRIB rsl0830963 (CIR 

P = 4.3 X 10"^^; DI P = 3.4 X 10"^^). AdditionaUy, we con- 
firmed that genetic variants in CDKALl (CIRP = 1.6 X 10"'^; 
DI P = 2.5 X 10"^^), HHEX (CIR P = 1.4 X 10"^; DI P = 
0.0021), GIPR (CIRP = 3.1 X 10"^; DIP = 0.0061), GCK(CIR 
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P = 0.0010; DI P = 9.2 X 10"^), KCNQl (CIR P = 0.0032; DI 
P = 0.0054), CDKN2A/2B (CIR P = 0.016; DIP = 0.0072), 
i^AD<Si (CIR P = 6.6 X 10"^), IGF2BP2 (DI P = 0.0013), 
JAZFl (DI P = 0.015), KCNJll (CIR P = 0.019), PPARG (DI 
P = 0.023), ZFAND6 (CIR P = 0.026), CHCHD2P9 (DI P = 
0.039), ARAPl (CIR P = 0.039), and BCLllA (CIR P = 0.048) 
were associated with reduced insulin secretion (Supple- 
mentary Table 2). 

Insulin action. The variants in MTNRIB (P = 6.3 X 10"^), 
CHCHD2P9 (P = 0.011), FADSl (P = 0.013), PPARG (P = 
0.041), and CDKALl (P = 0.042) were associated with re- 
duced insulin sensitivity (Supplementary Table 3). 
Glucagon secretion. To explore whether these genetic 
variants would also influence a-cell function, we analyzed 
glucagon concentrations during OGTT. The variants in 
BCLllA (P = 0.0069), HHEX (P = 0.0096), ZBED3 (P = 
0.033), and HNFIA (P = 0.047) were associated with ele- 
vated fasting glucagon levels, whereas variants in CRY2 
(P = 0.0021), IGF2BP2 (P = 0.0036), and TSPAN8 (P = 
0.034) were associated with reduced fasting glucagon 
levels. Furthermore, variants in IGFl (P = 0.0048), ZBED3 
(P = 0.0094), and N0TCH2 (P = 0.016) were associated 
with increased postprandial glucagon levels and the vari- 
ant in KCNJll (P = 0.027) with decreased postprandial 
glucagon levels (Supplementary Table 4). 
Effect of SNPs on relationships between metabolic 
variables 

Insulin secretion vs. sensitivity (CIR vs. ISI). We 

plotted the effects of the tested genetic variants on CIR 
and ISI (Fig. 1). Most SNPs were related to insulin secre- 
tion and sensitivity in a linear fashion, with SNPs in the 
MTNRIB, CDKALl, GCK, CDKN2A/2B, and PPARG loci 
in the lower left quadrant, reflecting both impaired insulin 
secretion and action, with MTNRIB being the most ex- 
treme. In contrast, variants in GIPR, FADSl, CAMKID, and 
TCF7L2 showed relatively high insulin sensitivity despite 
low insulin secretion. A genotype-centric analysis (Supple- 
mentary Fig. 1) also showed that homozygous carriers of 
MTNRIB, CDKALl, and GCK risk genotypes had CIR val- 
ues below the values expected from ISI values, whereas TT 
carriers of TCF7L2 showed rather high insulin sensitivity. 



Glucose vs. insulin levels. The variant in MTNRIB 
showed both high fasting glucose and insulin levels (Fig. 2 
and Supplementary Tables 5 and 6). At the same glucose 
level, CHCHD2P9 showed the highest and FADSl the 
lowest fasting insulin concentration, whereas at the same 
fasting insulin concentration, FTO showed the lowest and 
GCK the highest fasting glucose levels. Homozygous car- 
riers of the TCF7L2 risk variant had the lowest insulin 
levels, whereas homozygous carriers of HMGA2 showed 
the highest insulin levels (Supplementary Fig. 2). These 
differences were also partially seen at 2 h of the OGTT. 
Glucose vs. glucagon levels. The highest fasting glucose 
level related to glucagon level was seen for MTNRIB (Fig. 
3 and Supplementary Fig. 3), which was also reflected by 
high fasting and 2-h insulin-to-glucose ratios (Supplemen- 
tary Table 7 and Supplementary Figs. 4 and 5). In line with 
previous reports, the TCF7L2 variant showed low gluca- 
gon levels (9) and homozygous T-allele carriers showed 
high fasting and 2-h glucose levels, despite low insulin-to- 
glucagon ratio. 

Glucagon vs. insulin levels in human pancreatic 
islets. The risk allele carriers of variants in CAMKID, 
TCF7L2, GIPR, TSPAN8, and KLF14 showed lower glu- 
cagon levels whereas the risk allele carriers of CDKALl, 
BCLllA, CRY2, FADSl, N0TCH2, HMGA2 and GCK 
showed higher glucagon levels than expected from the 95% 
CI of the relationship between mean values of insulin and 
glucagon concentrations measured at normal glucose lev- 
els. At high glucose, carriers of the risk alleles at DGKB, 
GIPR, JAZFl, CDKN2A/2B, HMGA2, and N0TCH2 dis- 
played high glucagon levels, whereas those of TCF7L2, 
CAMKID, KLF14, FTO, and ZBED3 showed lower gluca- 
gon levels (Fig. 4). 

DISCUSSION 

In the current study, we demonstrate that the effects of 
some of the genetic variants associated with type 2 di- 
abetes and/or glycemic traits seem to involve the cross- 
talk between hormone secretion from the p- and a-cells. 
The novel findings in the current study included the ob- 
served associations of variants in 10 loci with glucagon 
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FIG. 1. Graphic representation of effects of genetic variants on the relationship between insulin secretion and insulin sensitivity in nondiabetic 
participants in the PPP-Botnia Study. Effect sizes are betas from linear regression adjusted for age, sex, and BMI. Whereas most genetic variants 
either increased insulin secretion or insulin sensitivity, those in MTNRIB, CDKALl, and GCK were located to the lower left quadrant, reflecting 
both impaired insulin secretion and insulin sensitivity, with MTNRIB being the most extreme. In contrast, GIPR and FADSl showed relatively high 
insulin sensitivity despite low insulin secretion. KCNQl t represents rs2237895, and KCNQl* represents rs231362. Black circles, P < 0.05; white 
circles, P > 0.05. 
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FIG. 2. Graphic representation of effects of genetic variants on the relationship between glucose and insulin levels during OGTT in nondiabetic 
participants of the PPP-Botnia Study. Effect sizes are betas from linear regression adjusted for age, sex, and BMI. KCNQIf represents rs2237895, 
and KCNQl* represents rs231362. Black circles, P < 0.05; white circles, P > 0.05. A: The relationship between fasting insulin and fasting glucose. 
MTNRIB showed the highest effect on fasting glucose despite elevated insulin levels. CHCHD2P9 showed increased whereas FADSl showed 
decreased fasting insulin. B: The relationship between 2-h insulin and 2-h glucose. GCK, CDKALl, ZBED3, and MTNRIB all showed high post- 
prandial glucose levels despite elevated insulin levels. 



concentrations, of which 6 loci (HHEX, IGF2BP2, CDKALl, 
FADSl, GCK, and CDKN2AJ2E) were also associated with 
reduced insulin secretion. 

Effects of variants on the relationship between CIR and 
ISI were similar to those observed by Grarup et al. (10), 
but we also showed that MTNRIB and CDKALl were 
clear outliers with strong effects on both impaired insulin 
secretion and sensitivity. These analyses used p-values 
that reflect fractions of the SDs. We therefore extended 
the analyses by analyzing the relationship between the 
absolute values of CIR and ISI in different genotype car- 
riers, which allowed us to obtain a quantitative effect of 
SNPs. These analyses differ from those presented by 



Grarup et al. (10) as we now related the effect of the ho- 
mozygous risk genotypes to the mean values of a reference 
population instead of expressing them as allelic changes in 
fraction of SDs. It can, of course, be questioned whether 
homozygosity for risk alleles would best describe their ef- 
fect on this relationship, but it has the advantage that it 
allows us to show the absolute effects of the variants. 

Importantly, we show that many of the genetic variants 
seem to affect the relationship between insulin/glucagon 
secretion and blood glucose levels both at fasting and 2 h 
during OGTT. This would imply that in the risk carriers, 
reduced early insulin secretion in conjunction with im- 
paired suppression of glucagon secretion and/or action 
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FIG. 3. Graphic representation of effects of genetic variants on the relationship between glucose and glucagon levels during OGTT in nondiabetic 
participants of the PPP-Botnia Study. Effect sizes are betas from linear regression a^usted for age, sex, and BMI. KCNQIf represents rs2237895, 
and KCNQl * represents rs231362. Black circles, P < 0.05; white circles, P > 0.05. ^: The relationship between fasting glucagon and fasting glucose. 
The highest fasting glucose level related to glucagon level was seen for MTNRIB. IGF2BP2 and CRY2 showed reduced glucagon secretion, whereas 
BCLllA and HHEX were associated with increased fasting glucagon. B: The relationship between 2-h glucagon and 2-h glucose. ZBED3, IGFl , and 
NOTCH2 showed elevated whereas KCNJll showed decreased postprandial glucagon secretion. 
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FIG. 4. Relationship between glucagon and insulin secretion in human pancreatic islets at low and high glucose levels. A: Secretion at low glucose 
(1 mmol/L). The risk allele carriers of variants in CAMKID, TCF7L2, GIPR, TSPAN8, and KLF14 showed lower glucagon secretion, whereas the 
risk allele carriers of CDKALl, BCLllA, CRY2, FADSl, NOTCH2, HMGA2, and GCK showed glucagon levels above the expected 95% CI of the 
relationship between insulin and glucagon secretion. B: Secretion at high glucose (16.7 mmol/L). The risk allele carriers otDGKB, GIPR, JAZFl, 
CDKN2AI2B, HMGA2, and NOTCH2 remained high glucagon secretion, whereas those of FTO, CAMKID, TCF7L2, KLFU, and ZBED3 had low 
glucagon. 
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would lead to increased risk of type 2 diabetes. An in- 
teresting finding was that, again, MTNRIB was different 
from the rest of the variants and showed elevated glucose 
despite high insulin-to-glucagon ratio. There is no obvious 
explanation for these findings apart from that the ratio 
between insulin and glucagon may be more important than 
the individual insulin and glucagon levels for the glucose- 
raising effect of MTNRIB. Also, MTNRIB is mostly 
expressed in (3-cells, whereas MTNRIA predominates in 
a-ceUs. Until now, no variant in the MTNRIA gene has been 
associated with type 2 diabetes or elevated glucose levels. 

In keeping with a previous study (11), the rs7903146 
variant in TCF7L2 was associated with lower glucagon 
values both in vivo and in vitro. In contrast to some pre- 
vious studies (12,13), we did not observe a significant re- 
duction in fasting or glucose-stimulated insulin secretion 
for this variant in the current study. The PPP-Botnia Study 
is a population-based study including people 18-75 years 
of age, most of whom had normal glucose tolerance, in the 
previous studies, the effect on insulin secretion was clearly 
stronger in hyperglycemic individuals or in those who 
subsequently developed diabetes. The same was seen in 
the current study; if we restrict the analysis to hyperglycemic 
individuals (fasting plasma glucose >5.5 mmol^), we see 
a clear effect of the TCF7L2 SNP on reduced insulin secre- 
tion (n = 1,440; (3 = -0.109 [0.035]; P = 0.002). This effect 
might be too small to fully explain the glucose-increasing 
effect of the SNP, nor is it easy to see an effect of reduced 
glucagon levels. One possibility is that the TCF7L2 variant 
has effects on hepatic glucose production as previously 
shown (13-15), which cannot be explained by effects of 
glucagon on the liver. 

From the figures relating the mean values of fasting and 
2-h insulin/glucagon with fasting and 2-h glucose con- 
centrations (Supplementary Fig. 6), it was obvious that 
homozygous carriers of HMGA2 differed most from the 
others. Thus, although carriers of two HMGA2 risk alleles 
showed increased insulin secretion and decreased insulin 
sensitivity, possibly protecting them from type 2 diabetes, 
they demonstrated low insulin/glucagon, which in turn led 
to elevated glucose. In vitro, homozygous risk allele car- 
riers showed a reduced suppressive effect of high glucose 
on glucagon secretion. HMGA2 encodes for a nonhistone 
protein with structural DNA-binding domains and may act 
transcriptional regulating factor. The current data 
suggest that it may involve a complex interplay between 
islet function and insulin action. 

The current study also provides some novel insight for 
the potential mechanism by which the BCLllA variant 
could increase the risk of type 2 diabetes. BCLllA rs243021 
was originally identified as a risk variant for type 2 diabetes 
in the DIAGRAM+ meta-analysis (1). The protein encoded 
by BCLllA is a zinc-finger protein with highest e3q)ression 
in the fetal brain and the adult central nervous system. Al- 
though its biological function is not fully elucidated, it might 
be involved in the signal transduction from cell to cell or in 
the neuronal regulation of hormonal secretion. Association 
of variants in BCLllA with reduced insulin secretion has 
previously been reported (16). Here we report novel asso- 
ciations of SL BCLllA variant with increased glucagon levels 
in vivo. Additionally, in vitro, homozygous risk allele car- 
riers of the BCLllA variant had glucagon levels above the 
expected 95% CI of the relationship between insulin and 
glucagon secretion in human pancreatic islets. Together, 
results from the present paper and previous reports suggest 
that this variant might predispose to increased risk of type 2 
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diabetes via reduced insulin and increased glucagon con- 
centrations. 

In conclusion, here we provide a comprehensive descrip- 
tion of how common genetic variants associated with type 2 
diabetes and/or glycemic traits influence insulin and gluca- 
gon secretion in vivo and in vitro in human pancreatic islets. 
The data clearly emphasize the need to consider effects on 
a-cell function when studying mechanisms by which these 
variants increase risk of type 2 diabetes. 
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